Piezomagnetism and Stress Induced Paramagnetic Meissner Effect in Mechanically 

Loaded High-T c Granular Superconductors 
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Two novel phenomena in a weakly coupled granular superconductor under an applied stress are predicted which 
are based on recently suggested piezophase effect (a macroscopic quantum analog of the piezoelectric effect) in 
mechanically loaded grain boundary Josephson junctions. Namely, we consider the existence of stress induced 
paramagnetic moment in zero applied magnetic field (piezomagnetism) and its influence on a low-field magneti- 
zation (leading to a mechanically induced paramagnetic Meissner effect). The conditions under which these two 
effects can be experimentally measured in high-T c granular superconductors are discussed. 



Despite the fact that granular superconductors have 
been actively studied (both experimentally and theoret- 
ically) for decades, they continue contributing to the va- 
riety of intriguing and peculiar phenomena (both fun- 
damental and important for potential applications) pro- 
viding at the same time a useful tool for testing new 
theoretical concepts [0. To give just a few recent exam- 
ples, it is sufficient to mention paramagnetic Meissner ef- 
fect [|| (PME) originated from a cooperative behavior 
of weak-links mediated orbital moments and found to be 
responsible for unusual aging effects [§ in high-T c granu- 
lar superconductors (HTGS). Among others are also re- 
cently introduced thermophase [QJ|] and piezophase ||] 
effects suggesting, respectively, a direct influence of a 
thermal gradient and an applied stress on phase differ- 
ence between the adjacent grains. Besides, using a model 
of random overdamped Josephson junction arrays, two 
dual time-parity violating effects in HTGS have been 
predicted 0,0. Namely, an appearance of magnetic 
field induced electric polarization along with the con- 
comitant change of the junction capacitance (magneto- 
electric effect pi] ) and existence of electric field induced 
magnetization (converse magnetoelectric effect [0) via 
a Dzyaloshinski-Moria type interaction mechanism. 

In this Letter we discuss a possibility of two other in- 
teresting effects expected to occur in a granular material 
under sufficient mechanical loading. Specifically, we pre- 
dict the existence of stress induced paramagnetic moment 
in zero applied magnetic field (piezomagnetism) and its 
influence on a low-field magnetization (leading to a me- 
chanically induced PME). 

The possibility to observe tangible piezoeffects in me- 
chanically loaded grain boundary Josephson junctions 
(GBJJs) is based on the following arguments. It is well 
known ||-fL4j that the grain boundaries (GBs) are the 
natural sources of weak links (or GBJJs) in granular 
superconductors. Under plastic deformation, GBs were 
found [ jl5| to move rather rapidly via the movement of 
the grain boundary dislocations (GBDs) comprising these 
GBs. As a matter of fact, using the so-called method 
of acoustic emission, the plastic flow of GBDs with the 
maximum rate of vq — 1mm/ s has been registered fllfl] 



in YBCO ceramics at T = 77K under the external load 
of a = 10 7 N/m 2 . Using the above evidence, in Ref.9 a 
piezophase response of a single GBJJ (created by GBDs 
strain field acting as an insulating barrier of thickness I 
and height U in a SIS-type junction with the Josephson 
energy J oc e~ iv ^) to an externally applied mechani- 
cal stress was considered. The resulting stress-strain and 
stress-current diagrams were found || to exhibit a quasi- 
periodic (Fraunhofer-like) behavior typical for Joseph- 
son junctions (JJs). To understand how piezoeffects can 
manifest themselves through GBJJs, let us invoke an 
analogy with the so-called thermophase effect suggested 
originally by Guttman et al. (as a quantum mechani- 
cal alternative for the conventional thermoelectric effect) 
to occur in a single J J and later applied to HTGS || . In 
essence, the thermophase effect assumes a direct coupling 
between an applied temperature drop AT and the re- 
sulting phase difference A<fi through a J J. When a rather 
small temperature gradient is applied to a J J, an entropy- 
carrying normal current /„ = L n AT (where L n is the 
thermoelectric coefficient) is generated through such a 
junction. To satisfy the constraint dictated by the Meiss- 
ner effect, the resulting supercurrent I s = 7 c sin[A0] 
(with I c — 2eJ/h being the Josephson critical current) 
develops a phase difference through a weak link. In other 
words, the temperature gradient stimulates a supercon- 
ducting phase gradient which in turn drives the reverse 
supercurrent. The normal current is locally canceled by 
a counterflow of supercurrent, so that the total current 
through the junction I = I n + I s = 0. As a result, super- 
current / c sin[A0] = — 7„ = —L n AT generates a nonzero 
phase difference via a transient Seebeck thermoelectric 
field leading to the linear thermophase effect jjj Acj) = 
- arcsin(LtpAT) ~ -L tp AT with L tp = L n /I C (T). 

By analogy, we can introduce a piezophase effect (as 
a quantum alternative for the conventional piezoelectric 
effect) through a J J |J. Indeed, a linear conventional 
piezoelectric effect relates induced polarization P n to an 
applied strain e as fl7| P n = d n e, where d n is the piezo- 
electric coefficient. The corresponding normal piezocur- 
rent density is j n — dP n /dt — d n e where e(er) is a rate 
of plastic deformation (under an applied stress a) which 
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depends on the number of GBDs of density p and a mean 
dislocation rate Vd as follows |Q e(er) = bpvd{cr) (where 
b is the absolute value of the appropriate Burgers vec- 
tor). In turn, Vd — vo(a/a m ) with a m being the so- 
called ultimate stress. To meet the requirements imposed 
by the Meissner effect, in response to the induced nor- 
mal piezocurrent, the corresponding Josephson supercur- 
rent of density j s = dP s /dt = j c sin[A(/)] should emerge 
within the contact. Here P s = —2enb is the Cooper 
pair's induced polarization with n — N/V the pair num- 
ber density, and j c = 2ebJ/HV is the critical current 
density. The neutrality conditions (j n + j s = and 
P n + P s = const) will lead then to the linear piezophase 



the network of superconducting grains is under the in- 
fluence of an applied frustrating magnetic field H, the 
total phase difference through the contact reads (where 



AH, a) 



— (fij ARij) ■ H - Ba ■ fij ■ 



(2) 



effect Acj) = — arcsin[d pp e(a) 



dpp e(cr) (with d p 



d n /jc being the piezophase coefficient) and the concomi- 
tant change of the pair number density under an applied 
strain, viz., An(e) = d pn e with d pn = d n /2eb. Given 
the markedly different scales of stress induced changes in 
defect-free thin films |l9f| and weak-links-ridden ceram- 
ics |Q, it should be possible to experimentally register 
the suggested here piezophase effects. 

To adequately describe magnetic properties of a granu- 
lar superconductor, we employ a model of random three- 
dimensional (3D) overdamped Josephson junction array 
which is based on the well known tunneling Hamilto- 
nian lEll- 
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It is well known |l],|lg[ll]|24| that the self-induced 
Josephson fields can in principle be quite pronounced 
for large-size junctions even in zero applied magnetic 
fields. So, to safely neglect the influence of these effects 
in a real material, the corresponding Josephson penetra- 
tion length Aj must be much larger than the junction 
(or grain) size. Specifically, this condition will be sat- 
isfied for short junctions with the size d <C Aj, where 
A,/ = \J ^o/Anfiojc \l with Xl being the grain London 
penetration depth and j c its Josephson critical current 
density. In particular, since in HTGS Xl — 150nm, the 
above criterium will be rather well met for d ~ lam 
and j c ~ 10 4 A/m 2 which are the typical parameters for 
HTGS ceramics Likewise, to ensure the uniformity 
of the applied stress a, we also assume that d -C A ff , 
where A CT is a characteristic length over which a is kept 
homogeneous. 

When the Josephson supercurrent R = I c sin <fiij cir- 
culates around a set of grains (that form a random area 
plaquette), it induces a random magnetic moment u s of 
the Josephson network ]2jJ] 



where {i} = fi is a 3D lattice vector, N is the number 
of grains (or weak links), J{rij) is the Josephson cou- 
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pling energy with r 



fj the separation between 



the grains; the gauge invariant phase difference is defined which results in the stress induced net magnetization 



= 4>°a - Aij, where 



with tfii being 



as 

the phase of the superconducting order parameter, and 
= |^ f? A(r) ■ dl is the frustration parameter with 

A(r) the electromagnetic vector potential which involves 
both external fields and possible self-field effects (see be- 
low); $o = h/2e is the quantum of flux. 

In the present paper, we consider a long-range inter- 
action between grains [^| Jl0|Jll| , |2^| (with Jfaj) = J) and 
model the true short-range behavior of a HTGS sample 
through the randomness in the position of the supercon- 
ducting grains in the array (see below). For simplicity, 
we shall ignore the role of Coulomb interaction effects as- 
suming that the grain's charging energy E c <C J (where 
E c = e 2 /2C, with C the capacitance of the junction). As 
we shall see, this condition is reasonably satisfied for the 
effects under discussion. 

According to the above-discussed scenario, under an 
applied stress the superconducting phase difference will 
acquire an additional contribution 5<pij{o') = —Bo ■ r\j, 
where B = d n eo/o m j c b with eo = bpvo being the max- 
imum deformation rate and the other parameters de- 
fined earlier. If, in addition to the external loading, 



M S (H,3) = — < u s >= I dfijdRijf(fij,Rij)fl s , (4) 



where V is a sample's volume and / the joint probability 
distribution function (see below). To capture the very 
essence of the superconducting piezomagnetic effect, in 
what follows we assume for simplicity that an unloaded 
sample does not possess any spontaneous magnetization 
at zero magnetic field (that is M s (0, 0) = 0) and that its 
Meissner response to a small applied field H is purely dia- 
magnetic (that is M S (H, 0) ~ —H). According to Eq.(2), 
this condition implies = 2-Kra for the initial phase dif- 
ference with m = 0, ±1, ±2, ... Incidentally, this is also 
a requirement for current conservation at zero tempera- 
ture || . 

In order to obtain an explicit expression for the piezo- 
magnetization, we consider a site positional disorder that 
allows for small random radial displacements. Namely, 
the sites in a 3D cubic lattice are assumed to move 
from their equilibrium positions according to the nor- 
malized (separable) distribution function f(rijRij) = 
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fr(rij)fn(Rij). As usual |8| JlO| , |ll| , it can be shown that 
the main qualitative results of this paper do not de- 
pend on the particular choice of the probability distri- 
bution function. For simplicity here we assume an expo- 
nential distribution law for the distance between grains, 
fr(f) = f(x)f(y)f(z) with flxj) = (l/d)e-*i/ d , and 
some short range distribution for the dependence of the 
center-of-mass probability fii{R) (around some constant 
value D). While the specific form of the latter dis- 
tribution is not important for the effects under discus- 
sion, it is worthwhile to mention that the former dis- 
tribution function f r (r) reflects a short-range character 
of the Josephson coupling in granular superconductors 
where J[fij) — Je~ K ' rij . For isotropic arrangement 
of identical grains, with spacing d between the centers of 
adjacent grains, we have k — (4, 4, 4) and thus d is of 
the order of an average grain size. 

Taking the applied stress along the x-axis, a = (a, 0, 0), 
normally to the applied magnetic field H = (0, 0, H), we 
get finally 



M s (H,a) = -M Q 



HtotjH, ct)/H 
+ H 2 ot (H,<j)/H 2 } 2 ' 



(5) 



for the induced transverse magnetization (along the z- 
axis), where H tot (H,a) = H — H*(a) is the total mag- 
netic field with H*(a) = (er/<7o)i?o being a stress- induced 
contribution. Here, M = I C SN/V with S = irdD 
being a projected area around the Josephson contact, 
H = $q/S, and cr = Vm{jc/ jd)(b/d) with j d = d n e 
and eo = bpvo being the maximum values of the dislo- 
cation current density and the plastic deformation rate, 
respectively. 

Fig.l presents the stress induced magnetization at dif- 
ferent applied magnetic fields, calculated according to 
Eq.(5). As is seen, in practically zero magnetic field the 
piezomagnetization is purely paramagnetic (solid line), 
exhibiting a strong nonlinear behavior. With increasing 
the stress, it first increases reaching a maximum, and 
then rather rapidly dies away. Under the influence of 
small applied magnetic fields (dotted and dashed lines), 
the piezomagnetism turns diamagnetic (for low exter- 
nal stress) with its peak shifting toward higher load- 
ing. At the same time, Fig. 2 shows changes of the 
initial stress-free diamagnetic magnetization (solid line) 
under an applied stress. As we see, already relatively 
small values of an applied stress render a low field Meiss- 
ner phase strongly paramagnetic (dotted and dashed 
lines) simultaneously shifting the peak toward higher 
magnetic fields. According to Eq.(5), the initially dia- 
magnetic Meissner effect turns paramagnetic as soon as 
the piezomagnetic contribution H*{a) exceeds an ap- 
plied magnetic field H . To see whether this can actu- 
ally happen in a real material, let us estimate the typ- 
ical values of the piezomagnetic field H* . By defini- 
tion, H*(a) = (a/a m )(j d /j c )(d/b)H where H = $ /S 



is a characteristic magnetic field, and a m is an ulti- 
mate stress field. Typically for HTGS ceramics 
S ~ 10 /im 2 , leading to H ~ 1G. To estimate the needed 
value of the dislocation current density jd, we turn to 
the available experimental data. According to Ref.14, 
a rather strong polarization under compressive pressure 
a/a m ~ 0.1 was observed in YBCO ceramic samples at 
T = 77 K yielding d n = 10 2 G/m 2 for the piezoelectric co- 
efficient. Usually ||-0@|O|, for GBJJs e - lO" 2 ^ 1 , 
and b ~ lOnm leading to jd — d n io ~ lA/m 2 for the 
maximum dislocation current density. Using the typi- 
cal values of the critical current density j c — l0 4 A/m 2 
and grain size d ~ 1/im, we arrive at the following es- 
timate of the piezomagnetic field H* ~ \Q~ 2 H$. Thus, 
the predicted stress induced paramagnetic Meissner ef- 
fect (PME) should be observable for applied magnetic 
fields H ~ 10~ 2 -£/o — 0.01G which correspond to the re- 
gion where the original PME was first registered 
In turn, the piezoelectric coefficient d n is related to an 
effective charge Q in the GBJJ as @ d n = {Q/S){d/b) 2 . 
Given the above-obtained estimates, we get a reasonable 
value of Q — 10 _13 G for the charge accumulated at a 
GBJJ. It is interesting to notice that the above values of 
the aplied stress a and the resulting effective charge Q 
correspond (via the so-called electroplastic effect [E6( ) to 
an equivalent applied electric field E — b 2 a/Q ~ l(rV/m 
at which rather pronounced electric-field induced effects 
in HTGS were either observed (like an increase of the 
critical current in YBCO ceramics |27|]) or predicted to 
occur (like a converse magnetoelectric effect [p~T[ ) . 

In conclusion, let us briefly discuss the contribution of 
the so-called striction effects |2^] (which usually accom- 
pany any stress related changes). According to Ref.28 
the Josephson projected area S was found to slightly de- 
crease under pressure thus leading to some increase of 
the characteristic field Hq = $o/S. In view of Eq.(5), it 
means that a smaller compression stress will be needed 
to actually reverse the sign of the induced magnetization 
M s . Furthermore, if an unloaded granular superconduc- 
tor already exhibits the PME, due to the orbital cur- 
rents induced spontaneous magnetization resulting from 
an initial phase difference 4>% — 27rr in Eq.(2) with frac- 
tional r (in particular, r = 1/2 corresponds to the so- 
called (2~|§] 7r-type state), then according to our predic- 
tions this effect will either be further enhanced by apply- 
ing a compression (with a > 0) or will disappear under a 
strong enough extension (with a < 0) able to compensate 
the pre-existing effect. Given a very distinctive nonlin- 
ear character of M S (H, a) (see Figs.l and 2), the above- 
estimated range of accessible parameters suggests quite 
an optimistic possibility to observe the predicted effects 
experimentally either in HTGS ceramics or in a specially 
prepared system of arrays of superconducting grains. Fi- 
nally, it is worth noting that a rather strong nonlinear 
response of the transport properties in HgBaCaCuO 
ceramics was observed pof under compressive pressure 
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with a/a m ~ 0.8. Specifically, the critical current at 
a = 9kbar was found to be three times higher its value 
at a = l.bkbar, clearly indicating a weak-links-mediated 
origin of the phenomenon (in the best defect-free thin 
films this ratio never exceeds a few percents p9[). 
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FIG. 1. The induced magnetization M a /Mo as a function of 
the reduced applied stress c/ero, according to Eq.(5) for differ- 
ent values of reduced applied magnetic field: H/Ho = 0.001 
(solid line), H/H = 0.01 (dotted line), and H/H = 0.1 
(dashed line). 
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FIG. 2. The magnetization M 3 /Mo as a function of the 
reduced applied magnetic field H/Hq, according to Eq.(5) for 
different values of reduced applied stress: <t/<to = (solid 
line), <t/<7o = 0.01 (dotted line), and cr/ao = 0.1 (dashed 
line). 
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